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ABSTRACT 


The purpose of this thesis is to provide the student 
of queuing systems with a vehicle through which a better 
understanding of the interrelationships between stochastic 


processes and queuing systems can be achieved. 
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I. INTRODUCTION 


A queuing system is comprised of action facilities 
(servers) and waiting facilities (queues or lines). The 
relationship between the components can be of many forms, 
each with their own characteristics. Consider a tool 
crib in a machine shop where the machinists are obliged 
to check out their ООС: or consider the telephone system, 
EsEwhuch callers request communication capacity. Both 
situations involve requests for action from limited resources, 
in these cases, servers or telephone circuits. In order to 
allocate resources properly, it is necessary to understand 
the interaction between requests for action and the actions 
themselves. Individually these can be modeled as stochastic 
processes, but in order to predict flow through the queuing 
s Tem, rt is necessary to understand the interrelation 
between the patterns of event occurrences. When demands 
occur at regular intervals and actions are of fixed dura- 
tion, it is relatively easy to predict flow; if there are 
more requests per unit of time than actions over the same 
time span, the queue will grow without bound unless the 
event relationship is changed. When variation is added to 
Me timing of events, the flow is not as obvious. For 
example, if the average number of requests for action 
equals the average number of actions, per unit of time, but 
both occur with variability, then it would be useful to be 


able to predict such quantities as expected delay and action 





БІСІТІТІУ 1ііе СІШЕ. Queuing models are used for this 
purpose. 

The purpose of this thesis is to provide the student 
of queuing systems with a vehicle through which a better 
understanding of the interrelationships between stochastic 
processes and queuing systems can be achieved. 

This is done through a combination of a computer simula- 


tion and a computer-driven graphical display. 
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TII. QUEUING MODEL SIMULATION 


This computer program simulates queuing systems in two 
Bages. First, a next-event type computer simulation creates 
a time history of the queuing system. Second, a dynamic 
graphical display depicts the occurrence of the queuing 
muocess in real time. 

Five basic, academically interesting, queuing models, 
each with parameter options, can be simulated.  Queuing 
models are generally described in terms of the distribution 
Of request occurrences (arrival distribution), the distribu- 
tion of action occurrences (service distribution), and the 
number of action facilities (servers). Specifically, the 
Emus between request occurrences are taken to be indepen- 
dent and identically distributed realizations from the 
arrival distribution. The duration of each action occur- 
rence is also taken to be an independent and identically 
distributed realizatlon from the service distribution. 
Variations on these assumptions occur for some models and 
are described in Section A. 

The abbreviation G/G/C describes a queuing system in 
which the arrival and service distributions are general 
(unspecified) and there are C servers. 

iiemirvye Options in this simulation are: 

(1) G/G/C 
(2) G/G/C with losses 
(3) G/G/C with feedback 





(4) G/G/C with finite source 

(5) Two simultaneous G/G/C queues for comparison. 
The specific choices of arrival distribution and service 
distribution which can be simulated are fully described 


in Section B. 


A. THE QUEUING MODELS 


1. WE 





Picture 1 


This is the simplest of the queuing models considered 
here. Each arrival waits in the queue until a server is 
available, enters service for a calculated period of time, 
then departs the systen. 

This model may be simulated with from one to six 


BEIVETS, 
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2. G/G/C with Losses 





Picture 2 


Consider making a long distance telephone call on 
Mother's Day. If all the trunk lines are busy, your request 
for service will be rejected. The G/G/C model with losses 
reflects this by rejecting arrivals who find the queue 
full. The queue limit used is specified as an input 
parameter. Note that rejected arrivals are lost and are 
not assumed to return later. The fraction of all arrivals 
who are served (not rejected) is recorded for subsequent 


display. 
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5. СОО ОКУТ Н ЕТП т е боптсе 





Picture 3 


Consider a closed-circuit educational television 
System at a school. The finite population consists of all 
the televisions in the system. When one fails, it enters 
the queue to be serviced. After repair, it re-enters the 
population pending the next failure. 

This situation is modeled by assuming that each 
action source, upon entering the idle population, is 
assigned a waiting time from the arrival distribution. 
When this time expires, the source requests service, 
joining the queue if necessary. As each action source 
terminates its service, it returns to the idle population. 
The user specifies the number of servers (one to six) and 


the size of the population. 
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4. G/G/C with Feedback to the Queue 
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Picture 4 


This is an extension of the G/G/C model, with from 
one to six servers, allowing customers leaving service to 
depart the system or to return to the queue. At each 
servlce completion, the customer leaves the system with 
probability P or returns to the queue with probability 1-P. 


The user specifies P as an input option. 
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5. Two Simultaneous G/G/C Queues 





Picture 5 


This option allows comparison of two standard 
G/G/C queues, each with up to three servers. All parameters 
of each system are specified separately, and, if desired, 
identical arrival streams may be utilized. This allows 
comparisons between different arrival and/or service distri- 
butions and associated parameters, and between a different 
number of servers. Also, the effect of the various queuing 


disciplines can be compared. 


B. ARRIVAL/SERVICE DISTRIBUTIONS AND QUEUE DISCIPLINES 
As well as options already described, the following 
arrival and service distributions are available to the 


user of the program. Any of the following probability 
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AS tributions can be used as an arrival or service 
distribution. 
|] DIS tributrons 


a. К - ERLANG (exponential when K-1) 





K-1 J 
F(t) = 1 - E LC) г 
J=0 ' 
" 
RATE À = K 
cy = 1 
YK 
Input parameters are: 
(1) K 
(2) Rate - in terms of mean number of events 


per minute. 


b. Hyperexponential 


E 2 


where: 


(ED тает 


(2) CV = 


| 


(3) À = Rate is specified and then 


Ay > ho chosen to satisfy-- 


p шс, 


T5 





2%Уу2-а 


ED Ucet cS Pho 

(1) Rate, specified in terms of mean number 
of events per minute. 

(2) Coefficient of variation (must be greater 
than 1.0) 


Degenerate (Fixed Rate) 


DIRA 
CA 
I r A 
Inbut parameter: 


(1) Rate À specified in terms of mean number 


of events per minute. 


Queuing Disciplines 


For comparison purposes, there are four queuing 


disciplines to choose from. These determine the priority 


by which members of the queue will enter service. 


FIFO - the first into queue is the first out. 


Over all disciplines which are independent of service times, 


this one minimizes the variance of delay. 


LIFO - the last into the queue is the first 


Some systems, by their nature, function in this manner. 


SSTF - the queue member with shortest service 


time is served first; this minimizes the mean delay over 
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all other disciplines, but the variance of delay is greater 
mem that of FIFO, 

d. RANDOM DRAW - applies to a system in which 
insufficient information is known to implement any other 
discipline. 

3. Length of the Simulation 

The period of simulation is specified in minutes 

Dy the user. All queues are empty and servers idle at the 


start of the simulation. 
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MEA THET METHOD 


A. THE COMPUTER SIMULATION 

The computer simulation is a modified next event type. 
There are two types of events, arrivals to the queue and a 
completion of service. With the first there will be an 
associated entry to service if there is an idle server, 
and with the last an associated entry to service if the 
queue is not empty. 

The simulation is begun with the generation of the first 
arrival time and its associated service time. The times are 
generated via the transformation of random numbers to stochas- 
tic variates of the type specified. At the occurrence of 
an arrival, the next arrival time and associated service 
time are generated. When a server is available, a waiting 
member of the queue enters service with completion time 
equal to the sum of his assigned service time and the time 
he departs the queue. Thus with the clock starting at zero, 
a Sie history of events is created which tells what happened 
and when. 

At the same time, the appropriate values are retained 
teem which statistics of the simulation can later be 


calculated. 


B. THE GRAPHICAL DISPLAY 
At the completion of the simulation phase, the time 


wa 


history is depicted dynamically on the cathode ray tube 
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ШШ еа time. Pictures 1 to 5 show the form of the display, 
with the man-like images in motion along the appropriate 
path, and the number in the queue box showing the queue's 
current length. 

The timing and motion of the display are controlled 
by the analog computer via the hybrid connections to the 
digital computer. The first connection is for the purpose 
of timing until the next event is to be initiated. In this 
case, the analog computer integrates from zero to a preset 
mene Value лп volts, at the rate ot one volt per second. 
The interrupt to the appropriate subroutine then occurs, 
and the analog computer is reset with the next time value, 
after which the integration is reinitiated. 

The second connection causes the dynamic motion to 
occur by continually checking events in progress and causing 
the position of the corresponding images to be incremented 
mepreapriately; this interrogation and incrementation occurs 
approximately 10 times per second depending on the settings 
of the analog computer. Movement increments are between 
Me half inch and two inches (depending on the path), 
@ausing a motion of between five inches and 20 inches per 


second. 
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IV. STATISTICS FROM THE SIMULATION 


STATISTICS 
o2 


3.99 OC 
‚95 BC 
‚84 L.C 


.24 DV 
TI Бу 


Ov „Б7 


. 38 . DOM ШІ 
.9г 03 „35 
82 Пі. EP 


24 PL EI 
ENTER 1. CONTINUE | 
2. PARAMETER LISTING 





At the end of the simulation period, the following 


statistics are calculated and displayed, as shown above. 


AT = time - average length of the queue 
BT = time - average number of busy servers 
LT = time - average number in the system = QT + BT 


QC = customer - average length of the queue 


BC = customer - average number of busy servers 


LC = customer - average number in the system = QC + BC 
DC = customer - average delay in the queue 

SC = customer - average service time 

WC = customer - average wait in the system = DC + SC 
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DM 


DO 


DX 


QM 
QO 


QL 


DV 


SV 
QV 


NC 
PS 


maximum delay experienced by any customer 
probability that a customer's delay is greater than 
zero. 

probability that a customer's delay is greater than 


XD. (The value of XD is an input option.) 


maximum queue length reached during the simulation 
probability that the queue length will be greater 
than zero. 

probability that queue length will be greater than LQ 


(the value of LQ is an input option). 


standard deviation of the delay in queue 
standard deviation of the service times 


standard deviation of the queue length 


total number of arrivals 


percentage not lost (losses model only) 
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V. PARAMETER INPUT 


5 RUNTIME 
3 ІС! 
е 
1 
L 


IDIST 
KK { | 
RATE 38.88 38.88 14.88 
cv 1.88 ИГ 1.88 
IDSPLN 1 { 


LO 1 МАХО 
XD 1.88 NPOP 
IDUAL 1 P 


ENTER 1. CONTINUE 
2. STATISTICS 





Picture 7 


Parameters can be input by the card reader or by the 
teletype terminal. In either case, the namelist form is 
used, where the entry form is NAME = X, NAME being the 
variable title and X being the value being assigned. For 
real variables, the decimal point must be entered. Name- 


list input is terminated with a *. 


A. SAMPLE INPUT CARD 





IDEV = 2, ITYPE = 3, RATE(2) = 3.0 * 


бар 





B. NAMELIST INPUT VIA THE GRAPHICS TERMINAL 

If this mode of entry is specified, one variable at 
a time is entered, followed by a carriage return. Entries 
are displayed, as entered, on the graphics screen. Termina- 
tron occurs with the entry of a * followed by a carriage 


Return. 


C. INPUT VARIABLES 
The following conventions are employed for specifying 
the options and parameters to be used. 
HYPE = 1 gives G/G/C 
2 gives G/G/C with losses 
3 gives G/G/C with finite source 
4 gives G/G/C with feedback 
5 gives two simultaneous G/G/C queues 
RUNTIME = desired length of simulation (minutes) 
IC = number of servers for queue 1 
ШЕШІ = number of servers for queue 2 (only applicable if 


ITYPE = 5 is chosen) 


IDIST(I) = 1. K - ERLANG 
2. HYPEREXPONENTIAL 
3. DEGENERATE 
KK(I) = K in K-ERLANG (use K=1 for exponential) 
RATE (I) = Rate in number per minute 
О) = Coefficient of variation (specify for IDIST(I) = 2 


only. This variable is computed otherwise.) 
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where I = 1 for arrivals to queue 1 (except finite model) 
2 for arrivals to queue 2 or finite population 
idle parameters 
3 for service for queue 1 (except finite model) 
4 for service for queue 2 or finite population 


service parameters 


MUSPDLN(I) = 1. -FIFO 
ep Mh eG, 
3. Shortest service time first 
4. Random draw 

where I = 1 for queue 1 


2 for queue 2 


LQ parameter in the statistic P(Q > LQ) 


XD = parameter in the statistic P(wait > XD) 


IDUAL = 0. Independent arrival streams 
l. Identical arrival streams 


(applicable to ITYPE = 5 only) 


MAXQ = maximum queue size (for G/G/C with losses only) 


NPOP = population size (for finite source option only) 


P = completion probability (for feedback option only) 
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VI. ТНЕ PROGRAM 


The program requires the use of the XDS 9300 digital 
computer, the CI 5000 analog computer, and a graphics 
terminal using the GATED monitor. 

The analog computer wiring diagram necessary to drive 
the display sequence appears in APPENDIX A. Before 
beginning, set Delay Flop Switches DFOO, DFOl, and DFO2 
to .1 MS and set counter FOO to 1 and FO4 to 0. Then select 
RUN and DIGITAL COMPUTER on the analog control panel. This 
provides for movement interrupts to occur at a rate of 
100/11 hertz. If desired, the timing values will be dis- 
played on the RATIOMETER by selecting TRUNK 420 prior to 
DIGITAL COMPUTER. 

The program consists of a FORTRAN program and a META- 
SYMBOL 9300 program. Because the history is buffered out 
to the device assigned the logical unit value 7, an addi- 
опа control card is needed. 

Job cards oe deck run should be: 
JOB 
^AGT 
“ASSIGN 7=DFIA 


-FORTRAN LS,GO 
FORTRAN DECK 
META 9300 SI,LO,GO 


META-SYMBOL DECK 
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PROGRAM STEPS 
SET-UP ANALOG 
COMPUTER 
PLSETOGATIED 
GRAPHICS MONITOR] 
SPECTFY 
GRAPHICS TERMINAL 


INPUT PARAMETERS INPUT PARAMETERS 
VIA TERMINAL VIA CARDS 
DISELAY DISPLAY 
STATISTICS PARAMETERS 
DISPLAY 


ОАЕ ПЕС 


SIMULATION 

CREATE HISTORY 

QUEUING SYSTEM 

DYNAMIC DISPLAY 

WRITE PARAMETERS 

AND STATISTICS 
DISPLAY DISPLAY 
PARAMETERS STATISTICS 








APPENDIX A 


ANALOG SCHEMATIC 


SCOPE Y2 
10 CPS ms (б) 
Ф - N201 
DO B INT 
RESET au o 
© ENS | i 53 
1М5 
C11 
2 ОР E RESET, 
2 ® 
pu a 







100 CPS 


PSET O&A 


-] V/SE 
-100 с 
^ 10, 10, un 
1 1 С11 
DAC] 
E 
T420 
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